Impact of internal bremsstrahlung on the detection of 7-rays from neutralinos 
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We present a detailed study of the effect of internal bremsstrahlung photons in the context of 
the minimal supersymmetric standard models and their impact on 7-ray dark matter annihilation 
searches. We find that although this effect has to be included for the correct evaluation of fluxes 
of high energy photons from neutralino annihilation, its contribution is relevant only in models and 
at energies where the lines contribution is dominant over the secondary photons. Therefore, we 
find that the most optimistic supersymmetric scenarios for dark matter detection do not change 
significantly when including the internal bremsstrahlung. As an example, we review the 7-ray dark 
matter detection prospects of the Draco dwarf spheroidal galaxy for the MAGIC stereoscopic system 
and the CTA project. Though the fiux of high energy photons is enhanced by an order of magnitude 
in some regions of the parameter space, the expected fluxes are still much below the sensitivity of 
the instruments. 

PACS numbers: 95.35.+d, 95.55.Ka, 98.52.Wz, 12.60.Jv 

The minimal supersymmetric (SUSY) extension of the 
standard model (MSSM) provides a natural candidate 
for dark matter (DM) in the form of a neutral, stable 
Majorana fermion, the lightest neutralino. At present, 
large efforts are being carried out to detect this SUSY 
DM by different methods, see [l| for reviews. 

In the case of the current imaging atmospheric 
Cherenkov telescopes (lACTs), the searches are based 
on the detectability of 7-rays coming from the annihi- 
lation of the SUSY DM particles in the halo of galax- 
ies [2]. Neutralinos annihilate at the one loop level into 
photons through the processes Q XX ~^ 77 1 XX ~^ 
with almost monochromatic outgoing photons of ener- 
gies E~f ^ — rri^ / 4m^ , respectively. More- 
over, neutralino annihilation can produce a continuum 
spectrum of secondary photons from hadronization and 
decay of the annihilation products, mostly from neutral 
pion decay, which typically dominates over the number 
monochromatic 7's in a large portion of the parameter 
space. lACTs in operation like MAGIC, HESS, VER- 
ITAS [3| or satellites-based experiments like the Fermi 
satellite Q play a very important role in this kind of 
DM searches. For these experiments, dwarf spheroidal 
(dSph) galaxies around the Milky Way represent a good 
alternative target option to e.g. the Galactic Center, al- 
ready observed in 7-rays but with null DM detection so 
far [6|. Dsphs are DM dominated systems with inferred 
very high mass-to-light ratios, and most of them are ex- 
pected to be free from any other astrophysical source 
that might contribute to a possible 7-ray signal. There- 
fore, the detection of 7-rays from them would probably 
imply a successful DM annihilation detection. 

Some of the present authors calculated in Ref. 0] the 
expected 7-ray flux due to neutralino annihilation in the 
Draco dSph for a typical lACT above 100 GeV. Draco 
is located at 80 kpc and is one of the dwarfs with more 



observational constraints, which have helped to better 
determine its DM density profile. 

The MAGIC telescope has already observed Draco in 
7-rays in the context of DM searches Q, but found no 
gamma signal above an energy threshold of 140 GeV. As 
a consequence, an upper limit for the flux {2a level) was 
set to be 1.1 x lO^^^ph cm~^s~^, assuming a power-law 
with spectral index —1.5 and a point-like source. This 
upper limit is 0(10"^ — 10^) above the values predicted by 
those SUSY models used in their analysis and therefore 
no constraints could be put on the parameter space. Also 
the Fermi collaboration has recently reported their upper 
limits for a possible 7-ray annihilation signal from Draco 
at lower energies 0], given that no significant gamma 
emission was detected above 100 MeV. 

Recently, in Ref. [13], it was noted that the photons 
arising from internal bremsstrahlung (IB) in some re- 
gions of the parameter space can dominate the spec- 
trum at energies near the neutralino mass. IB [13, [nl IS 
commonly referred to the emission of additional photons 
in neutralino pair annihilation into charged final states, 
XX ^ XX J {X being a charged lepton, a quark, a W bo- 
son or a charged Higgs) which is an unavoidable electro- 
magnetic radiative correction. In the Feynman diagrams 
these photons can be attached to the external legs repre- 
senting final state charged particles or to the propagator 
of the virtual charged particle exchanged by neutralinos: 
the latter diagrams are at the origin of the hard photon 
spectrum of IB. 

The aim of this brief report is to quantitatively study 
the impact of the IB on the DM detection prospects for 
lACTs. As an example, we will revisit Draco updating 
the results obtained in Ref. f?\ , this time fully taking into 
account the IB contribution to the expected annihilation 
flux. We note that a study of Draco including IB was 
already done in Ref. p^ . but only for a few benchmark 
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points of the parameter space. Here, we will perform a 
wider exploration of the parameter space and will extract 
more general conclusions on the real importance of the 
IB for 7-ray DM searches. 

The expected flux of photons with energy above the 
threshold of the telescope, F{E^ > ii^th), is given by the 
product of the so-called astrophysical factor J(5') times 
the particle physics factor, namely $_pp(£'^ > E'th)- 
J{^) represents the integral of the square of the dark 
matter density pdm along the direction of observation ^ 
relative to the center of the DM halo, and depends on 
the PSF of the telescope. In the case of Draco, for in- 
stance, the authors in Ref. used a cusp and a core DM 
density profiles, built from the latest stellar kinematic ob- 
servations together with a rigorous method of removal of 
interloper stars. In the same work, they also stressed the 
important role of the PSF, which is directly related to 
the angular resolution of the lACT and becomes crucial 
for a correct interpretation of a possible gamma signal 
due to neutralino annihilation. However, for the sake of 
simplicity, we will use here the value of J(vl') integrated 
over the whole spatial extent of the source as the value 
of the astrophysical factor. This value, that does not de- 
pend on the PSF any longer, can be well approximated 
by J = ^J-jj-i Jy p^]\i{r) dV, with D the distance from 
the Earth to the center of the DM halo and r the galac- 
tocentric distance inside it. For Draco we take a value of 
J = 3.7 X 10^^ Gev^ cm~^, which was calculated using 
the cuspy DM density profile given in Ref. Q. 

As for ^pp{E^ > Eth), which in the following we 
call fsusY, it includes all the particle physics informa- 
tions and is made up by the contribution of the con- 
tinuum spectrum (secondaries and IB photons) and the 
monochromatic photons (lines): 

f susy — fcont f lines , 
J lines — ^ 2m J 2mJ ' 

Here dN^/dE.^ is the differential yield of photons per 
annihilation to the final state / with branching ratio 
Bf. The factor in parenthesis is thus n^{Ej > Eth), 
the total number of photons per annihilation with energy 
greater than the threshold energy, {a^^ v) is the thermal 
averaged total neutralino annihilation cross, (cr^-y v) and 
{(Jz-y "v) the cross sections for annihilation into lines and 
TO^ the neutralino mass. 

As in Ref. fj\ we assume that the neutralino is the 
main component of the DM present in the Universe 
with abundance inside the cosmologically favored interval 
0.09 < n^h"^ < 0.13 (the most recent WMAP ^ inter- 
val at 3cr is 0.094 < noMh"^ < 0.128). We further require 
that SUSY models satisfy the LEP bounds on Higgs and 
chargino masses, m/i > 114 GeV and > 103.5 GeV, 
and constraints from 6 — 37 as explained in Refs. 
We first consider minimal supergravity (mSUGRA) mod- 
els, where the soft terms of the MSSM are taken to be 



universal at the gauge unification scale Mqut- The effec- 
tive theory at energies below Mqjjt is thus determined 
by four universal parameters: the common scalar mass 
Too, the gauginos mass mi/ 2, the trilinear couplings Aq, 
and the ratio of the Higgs vacuum expectation values, 
tan /?. In addition, the minimization of the Higgs poten- 
tial leaves undetermined the sign of the Higgs mass pa- 
rameter /i that we take positive. For the numerical com- 
putation of IB effects we use DarkSusy 5.0.5 [l3| which 
relies on ISASUGRA 7. 78 for the renormalization group 
equation evolution of parameters. 

The amount of photons coming from IB depends on 
the annihilation channels and therefore it depends on 
the nature of the neutralino in every particular SUSY 
model. We thus display in Fig. [1] the differential yield of 
photons for some representative points of the mSUGRA 
parameter space specified in Table |T] where the distinct 
contributions to fsuSY integrating the number of pho- 
tons above Eth = 100 GeV can be read. Model A is in 
the stau coannihilation region of the mSUGRA param- 
eter space: the mass of the lightest stau is TOf = 195 
GeV very close to = 188 GeV. Neutralino pair anni- 
hilation in T'^r" mediated by i-channel exchange of stau 
has the highest annihilation cross section thus IB is rel- 
evant. Here fiB is the dominant contribution being 10 
and 4.4 times greater than fsec and funes- Model B is 
in the funnel or resonances region: the mass of the CP- 
odd neutral Higgs is tua = 1211 GeV while = 598 
GeV, thus niA — 2m^ and WMAP bounds are satisfied 
due to neutralino pair annihilation into fermion through 
s-channel exchange of heavy neutral Higgs bosons. In 
this case no photon line can be attached to the virtual 
particles and the IB yield is negligible. Model C is in the 
focus point or hyperbolic branch region. The mass of the 
lightest chargino is m^± = 212 GeV, not much bigger 
than = 163 GeV and neutralino pairs annihilate into 
W^W~ through i-channel chargino exchange. The IB 
yield is small because is not much greater than tovi/ 
and photons energy has a cut off which corresponds to 
the kinematic endpoint x = \ — mw'^ /rn'^ ~ 0.75. Here 
funes is bigger than fsec and //b. Model D is another 
example in the focus point region. The mass of the light- 
est chargino is m^± = 954 GeV, almost degenerate with 
TO^ = 918 GeV. Neutralino pairs annihilate into W^W" 
through t-channel chargino exchange as in C but in this 
case 3> mw thus IB photons contribution is more im- 
portant and have endpoint at the neutralino mass: here 
fsuSY is dominated by fsec- The point E is similar to 
point A but with even more pronounced IB contribution 
(here to/j — 113 GeV). fiB is a factor 50 and 5 greater 
than fsec and fimes respectively but the the total fsusY 
is the smaller one. 

These ratios are quite general as can be seen in Fig. [2] 
where we present a larger scan on the parameter space 
of fsuSY versus the threshold energy setting Aq — and 
two values of tan /3, 10 and 50. In panel (a) where points 
are in the stau coannihilation region we can appreciate 
the largest contribution of IB, as shown by the points A 
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TABLE I. mSUGRA models used in Fig. [T] The values of mo, mi/2, ^o, rn^i are in GeV, the sign of fi is positive. The units 
of (o"xx^) cm'^ s~^, those of the f's, defined in Eq. Q, are GeV"^ cm'^ s~^. 
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127 


459 
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0.1 


B 


52 
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1377 


725 
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0.72 
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10~^ 


0.72 
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17 


2200 


430 


805 


162.8 


0.098 


2225 


0.04 


0.06 


0.02 


0.12 
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51 


8940 


2218 


-4221 


918.2 


0.099 


1203 


0.3 


0.003 


0.017 


0.32 
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110 


530 


-600 


218.4 


0.1 


11.2 


0.0014 


0.014 


0.073 


0.088 




FIG. 1. Continuum 7-rays spectra for the models in Table [11 The solid lines indicate the total yield, the dot-dashed lines the 
contribution of secondary photons, while dashed lines correspond just to the IB contribution. 




FIG. 2. The particle physics factor fsuSY, Eq. |T|, versus Eth, energy threshold of the detector. We set Aq — and fj, positive 
while mo and mi/2 have values such that the mSUGRA point predicts the neutralino relic density inside the WMAP bounds 
satisfying all the phenomenological constraints. Also plotted are the predicted sensitivity lines of MAGIC II and CTA for Draco 
corresponding to 50 hours of observation time and a 5a detection level. 



and E. The absence of IB photons of the point B is ev- 
idenced by the panel (b) where points are in the funnel 
region, while the panels (c) and (d) have points mostly 
in the hyperbolic branch and share properties with the 
points C and D. We also plot in those panels the sensi- 



tivity lines for Draco of the MAGIC telescopes in stereo- 
scopic mode [1^ and of the CTA project [13 as given by 
Montecarlo simulations for 50 hours of observation time 
and a 5ct detection level. We see that even including IB 
for a typical value of the threshold Eth = 100 GeV, the 
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(GeV) 

FIG. 3. Scatter plot of n^{E-^ > 100 GeV)x {oxx^) versus for continuum photons in mSUGRA and general MSSM. Also 
plotted are the predicted sensitivity lines of MAGIC II and CTA for Draco assuming Eth ~ 100 GeV, 50 hours of observation 
time and a 5a detection level. 



theoretical predictions are at least three orders of mag- 
nitude below the detection limit. 

We further perform a scan on more general parameter 
space both in mSUGRA and in a general MSSM with 
random soft terms at the electroweak scale (see Ref. |2| 
for details on the scanned parameter space). In Fig. [3] 
we plot the quantity nj{Ej > Eth) x (""xx^) versus 
together with the MAGIC II and CTA sensitivity lines 
assuming Eth ~ 100 GeV: as above we see that a boost of 
at least three and two orders of magnitude respectively 
is needed to reach the detection line. We note, however, 
that the effect of substructures in the dwarf may enhance 
the annihilation flux importantly (see e.g. Ref. (Tsj). 

In summary, we have reported a detailed study of the 
internal bremsstrahlung contribution to the expected 7- 
ray flux from neutralino annihilation. We have found 
that although this effect has to be included for evalua- 
tion of fluxes of high energy photons from neutralino an- 
nihilation, its contribution is relevant only in models and 
at energies where the lines contribution is dominant over 
the secondary photons. As a result, the most optimistic 
particle physics scenarios for DM detection (which typi- 



cally correspond to those where most of the flux is given 
by secondary photons) will not change substantially. On 
the other hand, being typically the IB yield at most an 
order of magnitude greater than the lines yield, the net 
increase on absolute flux is of the same order. As an 
example of the impact of the IB on previous works on 
DM search, we recalculated the DM detection prospects 
of the Draco dwarf galaxy for the MAGIC telescopes, up- 
dating a previous work by some of the authors. We flnd 
that though the effect can rise the flux by e.g. an order 
of magnitude at 100 GeV, the predicted fluxes are still 
at least three orders of magnitude below the sensitivity 
of the instrument both in mSUGRA SUSY scenario and 
in the general MSSM. The same is applicable to other 
lACTs, given their roughly similar sensitivities at those 
energies. 
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